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Abstract 

Knowledge of in situ copepod diet diversity is crucial for accurately describing pelagic food web structure but is challenging 
to achieve due to lack of an easily applicable methodology. To enable analysis with whole copepod-derived DNAs, we 
developed a copepod-excluding 18S rDNA-based PGR protocol. Although it is effective in depressing amplification of 
copepod 18S rDNA, its applicability to detect diverse eukaryotes in both mono- and mixed-species has not been 
demonstrated. Besides, the protocol suffers from the problem that sequences from symbiotic ciliates are overrepresented in 
the retrieved 18S rDNA libraries. In this study, we designed a blocking primer to make a combined primer set (copepod/ 
symbiotic ciliate-excluding eukaryote-common: CEEC) to depress PGR amplification of symbiotic ciliate sequences while 
maximizing the range of eukaryotes amplified. We firstly examined the specificity and efficacy of GEEG by PGR-amplifying 
DNAs from 16 copepod species, 37 representative organisms that are potential prey of copepods and a natural 
microplankton sample, and then evaluated the efficiency in reconstructing diet composition by detecting the food of both 
lab-reared and field-collected copepods. Our results showed that the GEEG primer set can successfully amplify IBS rDNA 
from a wide range of isolated species and mixed-species samples while depressing amplification of that from copepod and 
targeted symbiotic ciliate, indicating the universality of GEEG in specifically detecting prey of copepods. All the 
predetermined food offered to copepods in the laboratory were successfully retrieved, suggesting that the GEEG-based 
protocol can accurately reconstruct the diets of copepods without interference of copepods and their associated ciliates 
present in the DNA samples. Our initial application to analyzing the food composition of field-collected copepods 
uncovered diverse prey species, including those currently known, and those that are unsuspected, as copepod prey. While 
testing is required, this protocol provides a useful strategy for depicting in situ dietary composition of copepods. 
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Introduction 

As tlie most numerous animals in marine ecosystem, copepods 
are critical link of primary production to higher trophic levels, and 
important driver of the marine biological pump [1]. Copepods can 
be herbivores, carnivores and omnivores, or can switch their 
trophic mode according to the relative abundances of the prey 
species. They can feed on a variety of prey belonging to diverse 
taxa and size categories, including phytoplankton, protozoans, 
eggs and larva of aquatic organisms, including those of copepods, 
and detritus [2-4]. Although copepods demonstrate remarkable 
versatility in their prey, they also exhibit specific feeding 



preferences among different prey species based on the traits of 
prey, such as motility, cell size and shape, nutritional value, 
dissolved chemical cues, and cell surface properties [5-8]. For 
example, both laboratory and field incubation studies have shown 
that copepods preferentially graze on ciliates and dinoflageUates 
when diverse foods are offered because they have higher 
nutritional quality than other prey species and are rich in 
polyunsaturated fatty acids (PUFA) and eicosapentaenoic acid 
(EPA) that influence the growth, survival and fecundity of 
copepods [9-10]. In addition, copepods can also discriminate 
between individuals of the same species with different properties. 
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Table 1. Information on collection of field samples for different purposes. 





station* 


Sampling date 


Type 


Pre-treatment 


Purpose 


AV 


31-Oct-2007 


Live copepods 


Gut evacuation and fixed 


Primer test 


DYW 


25-Apr-201 1 


Live copepods 


Gut evacuation and fixed 


Primer test 


DYW 


28-Jul-2012 


Live copepods 


Gut evacuation and fixed 


Primer test 


SYB 


20-Apr-2010 


Live copepods 


Gut evacuation and fixed 


Primer test 


PRE 


28-Aug-201 1 


Fixed water sample 




Primer test 


DYW 


28-Jul-2012 


Live copepods 


Gut evacuation 


Feeding experiment 


PRE 


28-Aug-201 1 


Fixed copepods 




Diets analysis 



AV, Avery Point {4ri8.917'N, 72=3.81 'W), Connecticut, USA; DYB, Daya Bay (22 36.274'N, 114"34.0'E), South China Sea, China; SYB, Sanya Bay {18"12.794'N, 
114'34.0'E), South China Sea, China; PRE, Pearl river estuary (22 7.022'N, 1 13 52.1 75'E), South China Sea, China. 
doi:l 0.1 371 /journal.pone.Ol 03528.t001 



including biochemical composition, growtli stage and nutritional 
quality [11]. 

Current knowledge of copepod feeding is largely derived from 
incubation experiments, which do not necessarily provide infor- 
mation on true diet composition of copepods at sea. Although 
natural dietary information can be obtained from gut content 
analysis of wild-caught copepods, the currendy available micro- 
scopic technique is not only time-consuming but also challenging 
because partially digested prey can be extremely difficult to 
identify. Pigment analysis has also been used to investigate prey 
diversity, but it is limited to phytoplankton and has low taxon- 
resolving power [12]. Stable isotope analysis of organic materials 
and fatty acid analysis are helpful in tracing sources of carbon or 
nitrogen and can provide information of diet categories or trophic 
level, but still could not reconstruct species composition of the diet 
[13-14]. 

Molecular techniques have been widely used for detecting the 
prey composition of aquatic predators in the natural environment 
due to its sensitivity, specificity and rapidness [15], among which 
polymerase chain reaction (PCR) assay has been successfully used 
to detect the food composition of marine invertebrates, such as 
crustacean (e.g. amphipods, krUl, copepods) and moUusca (e.g. 
bivalves) [16-18]. Although PCR assay has also been explored in 
copepod grazing research, most of these studies so far have focused 
on detecting predetermined prey species ingested by copepods 
[4,18-21] or targeted only one type of prey at a time [22], which 
stiU cannot provide a whole picture of in situ dietary composition 
of copepods. 

18S ribosomal RNA gene (18S rDNA) is widely used in PCR as 
a highly sensitive gene marker because it consists of multiple copies 
in the genomes of eukaryotic organisms [23]. Universal 18S rDNA 
primers (e.g. 18ScomFl/ 18ScomRl) have been proven useful in 
amplifying nearly all of the eukaryotes [24]; however, their 
application to analyses of food composition in a predator is 
constrained by the inevitable concurrent amplification of predator 
18S rDNA. In our own experiments using DNA from whole 
copepods (without dissecting guts) as template and 18ScomFl/ 
18ScomRl primer set for PCR amplification, the clone libraries 
were typically overwhelmed by 18S rDNA of the copepods. We 
recently designed a pair of copepod-excluding universal 18S 
rDNA primers (18S Non-copepod F2R2), and our initial test 
showed effective depression of copepod 18S rDNA and revealed 
several symbiotic cHiates associated with various species of 
copepods [25]. However, the efficiency of the primers to detect 
a broad range of eukaryotes has not been demonstrated. Besides, 
the consistent amplification of 18S rDNA from the diverse 
symbiotic cUiates associated with copepods casts some problems 



about the applicability of this primer set aimed at recovering all 
diet species of copepods. In this study, we designed a blocking 
primer which excludes sequences from those symbiotic cHiates, 
and combined it with 18S Non-copepod F2R2 to yield a copepod/ 
symbiotic cUiate-excluding, eukaryote-common primer set (18S 
Non-copepod F2R2+ blocking primer, named as CEEC). We 
firstly examined the specificity and efficacy of the CEEC primer 
set by testing 16 copepod species, 30 algal species from 6 phyla, 7 
other species from 4 eukaryotic lineages, and a field-collected 
microplankton sample, and then tested the applicability of the 
primer set to reconstructing diet composition for both lab-reared 
and field-coUected copepods. 

Materials and Methods 

Algal cultures 

Algal cultures were either isolated from South China Sea or 
obtained from the National Center of Marine Algae and Microbes 
(NCMA; formerly Provasoli-Guillard Center for the Culture of 
Marine Phytoplankton or CCMP), Bigelow Laboratory for Ocean 
Science, Maine, USA. Algae were cultured in f/2 or f/2-Si 
medium (pH 8.0, salinity 33± l%o) at 20°C on a 14:10 h light:dark 
cycle at ~100 |J.mol photons -m ^-s ' (Table SI). 

For each species, about 10,000 cells were harvested during the 
exponential growth phase by gentle filtration (<50 mm Hg) onto 
0.45/3-(im polycarbonate membrane, and then immersed in 
0.5 mL DNA lysis buffer (10 mM Tris-HCl pH 8.0, 0.5% SDS, 
100 mM EDTA pH 8.0, 200 |xg mL"' proteinase K) for 
subsequent DNA extraction. 

Individual species of other eukaryotes 

Seven other eukaryotic species were obtained from field 
seawater/aquaculture farms of Daya Bay (22°34.76'N, 
114°31.52'E) or coral ecosystem in Sanya Bay (18°12.996'N, 
109°27.003'E) (Table SI). Individuals or tissues from these 
eukaryotes were thoroughly ground on liquid nitrogen and then 
suspended in 0.5 mL lysis buffer for DNA extraction. 

Field-collected natural microplankton samples 

For field microplankton samples, 500 mL seawater was 
collected from 1-2 m below the surface in the Pearl River Estuary 
(22°7.022'N, 113°52.175'E) using a Niskin bottle (Table 1) and 
were fixed on site in neutral Utermohl's solution [26] at 2% final 
concentration for later analysis of ambient phytoplankton com- 
munities, a method proven to keep DNA intact for several months 
at room temperature [27]. Upon arrival at the laboratory, 50 mL 
subsamples were taken and settled down for two days in dark, and 
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Phaeodactylum tricornutum KF733528 
Nitzschia sp. IVID1 
Chaetoceros gracilis 

Unclassified diatom KF733539 

• Skeletonema costatum 

Thalassiosira weissflogii KF733529 
Nannochloropsis oceanica 
-Poteiioochmmonas sp. KF733553 

Gracilariopsis chorda KF733542 



Baelllariophyceae 



Chrysophyta 

Rhodophyta 



• Porphyridium cruentum 

Chlamydomonas sp. KF733532 
^^^y Platymonas helgolandica KF733530 

Tetraselmis subcordiformis KF733531 
_j- Unclassified chlomphyceae MC1 KF733533 CWoropliyta 

Scenedesmus sp. KF733545 
\t- Dunaliella sallna KF733535 

~ Haematococcus sp. KF733538 

^ Dysteiiidae sp. KJ569308 

• Euplotes rarlseta 



Icillophora 



■ Scrlppsiella trochoidea KF733540 
— Prorocentnim rhathymum KF733536 

Symbiodinium goreaui KF733550 

Prorocentnim lima KF733552 

- Alexandrium tamarense KF733551 



Coolia sp. KF733525 



-Ostreopsis sp. KF733537 



v-Amphidinium carterae KF733534 
' L Amphidinium sp. KF733526 
^ — Amphidinium klebsi KF733541 



I 
n 



^^^^^^^^ilsochrysis zhanjiangensis KF733547 

^ Pleurochrysis sp. KF733549 I Haptophyta 

^ • Isochrysis galbana \ 

Tubastraea coccinea KF733548 I Cnidaria 

Unclassified sagittidae KF733546 



Tegillarca granosa KF733527 
' Geloina erosa KF733543 

- Asaphis violascens KF733544 
■ Rat sp. M11 188.1 



I Chaetognatha 



Mollusca 



0.05 



>75 

50-75 
< 50 



Figure 1. Maximum Likelihood (IVIL) phylogenetic tree of 18S rDNA gene from different eukaryotes amplified by CEEC primer set. 

Sequences marked with solid circles are 100% identical to those under GenBank accession numbers AY229897.1, HM805045.1, JQ315726.1, 
HIVI149540.1, GQ246179.1, HM246242.1 and AJ305248.1. Other sequences obtained for mono-species were all submitted to GenBank and the 
accession numbers were shown in the tree (KF733525-KF733553, KJ569308). 
doi:1 0.1 371/journal.pone.01 03528.g001 



then concentrated to 1 niL which were then identified and 
enumerated in a Sedgwick-Rafter counting chamber under an 
Olympus BX51 microscope. Another 200 mL subsamples were 
taken and centrifuged at 3000 xg for 15 min, and the concentrated 
samples were pelleted by centrifugation at 12000xg for 5 min. 
The cell pellets were re-suspended in 0.5 mL DNA buffer for 
DNA extraction. 

Starved copepods 

To verify the non-amplification of 18S rDNA from copepods by 

CEEC primers, live copepods were collected using a zooplankton 
net (200-|J.m mesh) from the Avery Point campus of University of 
Connecticut in Long Island Sound, Pearl river estuary, Sanya Bay 
and Daya Bay in South China Sea (Table 1). No specific permits 
were required for the described field studies. A total of 16 species 
from 10 g(m(^ra w(^r(^ obtained (Table 2). They ^vl're kept alive in 
filtered si;awat('r during transportation to the laboratory and were 
starved in the laboratory to evacuate the gut contents, then 
identified and sorted under a stereoscope. The adult female 
individuals of the dominant species, Acartia erythraea, were sorted 
out for later feeding experiment (next section). Other copepods 
were homogenized for DNA extraction as reported [25] . 

Copepods fed on artificially mixed diet 

To test whether the primer set could efficiendy amplify 
individual species in a mixed diet, a short-time laboratory feeding 



experiment was carried out. Star\'ed copepods {A. erythraea) were 
introduced into the feeding chamber at 3.75 ind-mL~ and fed 
with a mixture of different algal species with equal biovolume basis 
(~1.0xlO'' |xm"^), iadudxag Prorocentrum donghaiense, Thalassio- 
sira weissflogii, Phaeodactylum tricornutum and Tetraselmis 
suecica. Every 3 minutes, copepod samples were taken and 
checked under a stereoscope to determine the feeding condition. 
Once the guts of copepods were filled and no fixes ^vere detected 
(< 15 min), samples were collected gentiy and rinsed carefully with 
filtered seawater, and then preserved for DNA extraction. 

Field-collected copepod samples 

Field copepod samples were also collected from the same 
sampling site in Pearl River Estuary waters (22°7.022'N, 
113°52.175'E) using a 505-|xm zooplankton net (Table 1). Once 
retrieved from the cod end, the sample was preserved immediately 
using neutral Utermohl's solution and transferred to laboratory 
where the dominated species, Canthocalanus pauper, was sorted 
out for DNA extraction. 

DNA extraction 

Algal and microzooplankton cell pellets and copepod homog- 
enates were incubated in lysis buffer for three days at 55''C for 
thorough cell lysis, and DNA was extracted following a modified 
CTAB protocol, and finally eluted in 50 ^lL 10 mM Tris-HCl 
(pH 8.0), as previously reported [24]. 
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Figure 2. 18S rDNA phylogram for grazer Canthocalanus pauper ambient planlcton community in Pear River Estuary. Maximum 
Lil<elihood (ML) tree was shown here and only representative clones from each major lineage were included in the tree. The color circles denote the 
sources of the clones. XZZ-W, ambient water sample; XZZ-F, Can. pauper sample. 
doi:l 0.1 371/journal.pone.01 03528.g002 



PCR protocol 

Due to the previously reported prevalent association of 
symbiotic ciliate with copepods [25], we developed a blocking 
primer (ciliatel8Sblk2: 5'-CCCAATCCTGATCCAGGGAGGT- 
AGTGACA.'\GAAATAACAACCTAGTCGCAAGACACACC- 
3') to specifically depress the amplification of 18S rDNA from 
those ciliates. Combined primers Non-copepod 18SF2, cilia- 
te 18Sblk2 and Non-copepod 18SR2 were used in a ratio of 
1:4:1. DNA from 37 eukaryotic organism, 16 copepod species, 
and a natural microplankton sample were PCR-amplified using 
this primer set under the same program as previously reported 
[25]. PCR products were purified and cloned into T-vector, and 
50 to > 100 clones were sequenced as reported [24]. To exclude 
the possibility that the negative result of the copepods DNAs 
was due to poor quality of the extracted DNA, each DNA 
sample was also PCR-amplified using the universal 18S rDNA 
primer set under the PCR conditions as reported previously 
[24]. 

Sequence data analysis and phylogenetic inference 

Sequences obtained were searched against the GenBank 
database using the Basic Local Alignment Search Tool. The best 
hits were collected and ahgned with the new sequences obtained in 
this study using CLUSTAL W (1.8) [28]. Maximum Likelihood 
(ML) tree was inferred from the ahgned dataset. 

Results 

Efficacy of CEEC primers to depress amplification of 
copepod 18S rDNA 

PCR with CEEC primer set on the DNA from 16 starved 
copepod species all gave negative results (Table 2). However, PCR 
with the universal 18S rDNA primer set all produced an amplicon 
with an expected product of ~ 1 .8 kb, indicating that the negative 
results with CEEC primers were not due to poor DNA quality, 
and cloning and sequencing results showed that they were all 
copepods or symbiotic ciliates. These results thus indicated that 
the primer set could eflFectively avoid amplification of 18S rDNA 
from the common copepod species and successfully block 
symbiotic ciliates sequences. 

PCR efficacy in detecting individual eukaryotic species 

The efficacy of CEEC primer set was confirmed by the positive 
PCR results from all 37 different species examined (Table SI), 
which represented a wide phylogenetic range of eukaryotes 
(Figure 1). Serjuences obtained for Chaeloceros gracilis, Skeleto- 
nema marinoi (coslalum), Nannochloris oculata, Isochrysis galbana 
and Nitzschia sp. (MDl) were 100% identical to the sequences 
from corresponding species documented in GenBank 
(AY229897.1, HM805045.1, JQ315726.1, HM149540.1 and 
GQ246179.1, respectively). Sequences from dinoflagellates Coolia 
sp., Amphidinium sj). and Prorocentrum sp., which were not 
morphologically identifu-d to species, gave 98-99% identity to 
sequences from the respecti\ e genera in GenBank. All other algal 
sequences showed 99% identity to those from different strains of 
the same genus or family. Sequences from ciliate Dysteriidae sp. 
and Euplole.s sp. showed 99% and 100% identity to a cyrtophorid 
ciliate Dysleriidae sp. (KF384515.1) and Euplotes rariseta 
(AJ305248. 1), respectively. Sequences from the other metazoan 



species all gave 98-99% identity to those from corresponding 
species documented in the database. 

PCR efficacy in revealing natural microplankton 
community 

Microscopic observation of water sample showed 28 different 
species, including 19 diatoms, 8 dinoflagellates and 1 metazoan 
species (Appendicularia: Oikopleura sp.). Of the dominant diatom 

(70%), Rhizosolenia showed the highest diversity, with 7 species in 
this genus, followed by Coscinodiscus, which contained 3 species. 
Gymnodiniaceae and Prorocentraceae dominated the dinoflageUate 
community (Table 3). Concerning the cell abundance, diatoms 
dominated the plankton community (98%), among which 
Rhizosolenia sp. was the most abundant species and accounted 
for about 51% of total abundance, followed by Thalassiosira sp. 
with 15%. 

Seventy six clones of 18S rDNA PCR product were sequenced 
for the water sample, of which 35 genotypes were resulted 
(Figure 2), and as BLAST analysis showed, these comprised 18 
species of diatoms, 9 species of dinoflagellates, 1 species of 
Cryptophyta, 1 species of Chlorophyta, 1 specic-s of Chrysophyta, 
1 species of Cladocera, 2 species of Appendicularia and 2 
unclassified nanoflagellates. Sequences from diatom dominated 
the clone library (47% by clone numbers), with Rhizosolenia sp. 
being the most abundant species. For the phytoplankton 
community, diatom species from Rhizosolenia, Thalassiosira, 
Coscinodiscus, Biddulphia, Skelelonema and Nilzschia, and dino- 
flagellates from Gymnodiniaceae and Prorocentraceae, were both 
detected in microscopic results and molecular detection. Com- 
pared to the microscopic results, a higher diversity of species was 
detected by the PCR assay (Table 3). These indicated that CEEC 
primers can successfully be used in amplification of DNA from 
assemblage of different organisms except those from copepods and 
symbiotic cUiates. 

PCR efficacy in reconstructing diet composition of 
experimentally fed copepods 

BLAST analysis of the sequences obtained from lab-fed 
copepods against the GenBank database showed that they were 
highly similar (>99% identity) to sequences from T. weissflogii 
(GU594641.1), P. Iricornulum (FR744760. 1), Pro. donghaiense 
(DQ336054.1) and T. suecica (FJ559381.1) in GenBank. This 
result indicated that the primer set could amplify all the mixed 
prey species in the guts of the laboratory fed copepods. 

Analysis of in situ copepod diets 

For the field-caught copepod (Can. pauper) samples, 34 
sequences were obtained, of which 6 taxa were identified 
(Figure 3), including those that were related to dinoflagellates, 
cryptophyte, radiolaria, echinoderm and tunicate (Table 4, 
Figure 2). Dinoflagellates were the dominant species in the 
copepod gut with one sequence showing high similarity to 
Karlodinium veneficum [[N986577, 99% identity) and most of 
other sequences relatively closely (90%-93% identity) related to 
documented Karenia mikimotoi (FR865627.1). The other 5 
sequences were close to that of the brittle star Ophionereis 
reticulata (DQP60805, 97% identity), the tunicate Oikopleura sp. 
(AB013015, 99% identity), the cryptophyte Goniomonas sp. 
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Table 3. Species composition of water sample by microscopic and molecular analyses. 



Genus 


Microscopic analysis 






Molecular analysis 




Cell density 
(cells/L) 


Percentage of 
abundance (%) 


Species number 


Species percentage (%) 


Clone number 


Percentage (%) 


Diatoms 


Rhizosolenia 


1.48x10'* 


56.68 


7 


25.93 


10 


13.17 


Thalassiosira 


3.88x10^ 


14.90 


2 


7.41 


7 


9.20 


Coscinodiscus 


2.4x10'' 


0.92 


3 


11.11 


3 


3.95 


Biddutphia 


1.2x10" 


0.46 


2 


7.41 


2 


2.63 


Skeletonema 


1.64x10^ 


6.30 


2 


7.41 


2 


2.63 


Nitzschia 


4.84x10^ 


18.59 


1 


3.70 


4 


5.26 


Navicula 


8.0x10^ 


0.31 


2 


7.41 






Unclassified 










8 


10.53 


Dinoflagellates 


Gymnodiniaceae 


1.6x10" 


0.61 


2 


7.41 


6 


7.90 


Prorocentraceae 


1.2x10'' 


0.46 


2 


7.40 


2 


2.63 


Ceratium 


8.0x10^ 


0.31 


1 


3.70 






Peridiniaceae 


8.0x10^ 


0.31 


2 


7.41 


1 


1.32 


Heterocapso 


4.0x10^ 


0.15 


1 


3.70 






Syndiniales 










7 


9.20 


Nanoflagellates 










2 


2.63 


Cryptophyta 










1 


1.32 


Chlorophyta 










1 


1.32 


Chrysophyta 










2 


2.63 


Metazoan 


35(ind/L) 








18 


23.68 


SUM 




100 






76 


100 



doi:1 0.1 371 /journal.pone.Ol 03528.t003 



(AY360454, 91% identity) and the radiolarian Acanthostamus 
purpurascens {JN811224, 99% identity), respectively. 

Discussion 

Molecular techniques have the advantages for reconstructing 
diet composition of natural population of copepods and other 
zooplankton in that they can detect partially digested, morpho- 
logically indistinguishable, and low-abundance species. Yet their 
applications have been restricted by the lack of broad-target DNA 
markers [17]. The present study was an attempt to develop and 
vahdate a copepods/symbiotic ciliates-excluding universal 18S 
rDNA-based PCR assay. We systematically evaluated the utility of 
this PCR assay for detecting diverse diets of copepods in natural 
environments by testing it against individual species of eukaryotic 
organisms, mixed species, natural assemblage, laboratory fed 
copepods, and field-caught copepods. The results indicate that this 
PCR assay largely meets the requirement of an easUy accessible, 
universal, and sensitive method. 

One potential major issue of applying such a PCR assay to 
detect diet diversity in copepods is the interference of copepod 
DNA, which is usually dominant in the DNA sample extracted 
from whole copepods [29]. Although isolating guts for DNA 
extraction can reduce this problem, it is time-consuming and still 
cannot guarantee complete removal of gut tissue-derived copepod 
DNA. In the present study, we showed that use of CEEC primers 
was very helpful in depressing amplification of copepod 18S rDNA 
for all copepods we examined. Most of the copepods tested were 



dominant species of coastal ecosystem from temperate or 
subtropical/ tropical regions, such as A. tonsa in coastal of Atlantic 
Ocean [30] and A. erythraea in South China Sea [31]. These 
copepods had different individual size, from<l |j,m (such as P. 
parvus) to>3 |im (such as S. suhcrassus) [31] and exhibited 
different feeding habits (herbivorous, carnivorous and omnivorous) 
[32]. Therefore, the CEEC primer set can be used to analyze diets 
for these, and likely other, copepods in various natural marine 
environment in future. The method will be particularly useful for 
small copepods, whose guts are very difficult to dissect for prey 
identification and feeding ecology is less understood. 

For a PCR assay to be useful in detecting food diversity, a wide 
taxonomic coverage is needed. In the present study, we used the 
CEEC primer set to amplify 37 different eukaiyotic species. Most 
of these were phytoplankton, which is considered as main diet 
components of mesozooplankton [22]. Because previous micro- 
scopic studies implied that the herbivores and omnivorous 
copepods might feed on unidentified diet materials derived from 
flagellates, aloricate cHiates, athecate, dinoflagellates or detritus 
[33], we also tested 2 species of free-living ciliates, which had been 
shown to be important prey item for copepods in laboratory 
experiments, and 5 field-collected marine metazoan, whose eggs, 
larvae or body parts might be ingested by copepods [3]. All the 
above phytoplankton and other eukaryotic species were success- 
fully amplified, which were verified by sequencing results, 
indicating that the PCR assay is efficient in detecting a wide 
range of different mono eukaryotic species. In addition, a natural 
microplankton sample from the ambient water was chosen to 
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Figure 3. Rarefaction curves for water sample and Canthocalanus pauper done libraries. The asymptote reached in XZZ-F suggests that the 
numbers of unique taxa in diets of Can. pauper was about 6, and no plateau reached for the XZZ-W plot indicates potentially higher diversity in water 
sample. 

doi:10.1371/journal.pone.0103528.g003 



evaluate whether the primer set can detect diverse species from a 
mixed sample. As the clone and sequencing results showed, much 
higher diversity of eukaryotic species (more than 35 species) was 
detected than that revealed by microscopic identification (28 
species). Almost all the phytoplankton groups revealed by 
microscopic identification were present in the clone library except 
Ceratium and Heterocapsa, which were missed likely due to 
insufficient sequencing depth, as suggested by the still growing 
taxon number in the rarefaction curve (Figure 3). It is noteworthy 
that nanoflagellates and cryptophyta, which might be diflicult to 
identify under light microscope, were detected, indicating its 
higher sensitivity to allow detection of small-sized and low- 
abundance taxa [30] . Taken together, all our results show that the 
primer set will be useful in reveahng in situ copepod diet diversity. 

The next challenge is to verify that the PCR assay would be 
sensitive and specific enough to detect prey in the gut when DNA 
from whole copepods is used. When our PCR protocol was 
applied to detect prey species in laboratory-fed copepods, all the 
prey species offered were successfully amplified. In the application 
to field-collected copepods, diverse and some unsuspected species 
in the diets were detected. Our results showed that the diet of Can. 



pauper' ranged from phytoplankton to metazoan, which is 
consistent with the previous report that Can. pauper is omnivorous 
[2] . It is reported that dinofalgeUates contribute a large portion of 
the carbon in the diets of common estuarine species [34]. In 
accordance, with the protocol developed here, we found many 
dinoflageUate 18S rDNA sequences. As some species of crypto- 
phyte (such as Rhodomonas sp.) has been considered to be prey of 
copepods, we also detected cryptophyte although the species 
(Goniomonas sp.) has not been reported [19]. Skeletal remnants of 
radiolarians have also been found in the guts of large copepods 
(e.g. Pleuromamma xiphias and Euchirella messinensis) by gut 
contents analysis in Sargasso Sea [3]. As gut dissection for small 
copepods, such as Can. pauper studied here, is difficult, metazoans 
as prey for copepods were rarely reported. Furthermore, brittle 
star and tunicate have not been reported as prey of copepods, but 
retrieved from field-caught copepods in this study. Therefore, the 
PCR protocol developed in this study has higher sensitivity and 
resolution in detecting in situ food composition of copepods. 

Besides excluding copepods, our PCR protocol was also 
designed to exclude epibiotic apostomatid cHiates (mainly 
Vampyrophrya-like and Vorticella-like species) by using a blocking 



Table 4. Taxonomic distribution of 18S rDNA clones retrieved from Can. Pauper. 



Clone ID 


Best hit species 


Best hit Accession No. 


E-value 


identities 


Clone number 


Class 


XZZ-F-1 


Karenia mikimoto 


FR865627.1 


0 


90%-93% 


29 


Dinophyceae 


XZZ-F-4 


Ophionereis reticulata 


DQ060805 


0 


97% 




Echinoderm 


XZZ-F-5 


Oikopleura sp. 


AB013015 


0 


99% 




Appendicularia 


XZZ-F-6 


Goniomonas sp. ATCC 50108 


AY360454.1 


0 


91% 




Cryptophyte 


XZZ-F-7 


Acanthostaurus purpurascens 


JN8n224 


0 


99% 




Radiolaria 


XZZ-F-8 


Karlodinium veneficum 


JN986577 


0 


99% 




Dinophyceae 


SUM 










34 





doi:l 0.1 371 /journal.pone.Ol 03528.t004 
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primer. However, it is important to verify that the protocol would 
not block out other ciliates that might be food of copepods. Two 
free-living ciliates, which are considered to be important prey item 
for copepods [2], could be successfully amplified by the primer set. 
Similarly, even though the protocol was intended to exclude grazer 
copepods, we did not want to also excludes othc-r crustaceans. From 
the ambient water sample, we detected Cladocera (95-96% 
identical to Cercopagis pengoi, EF189620.1), which accounted for 
17% of total clones. Furthermore, some other species (such as 
Farfantepenaeus duorarum) were also detected from copepod 
samples (unpublished data), indicating that the primer set could 
specifically evade amplification of copepod but other crustacean, 
larvae or body parts of which might be potential preys of some 
copepods [3]. 

Compared with pigment analysis or gut contents microscopy, 
this method can reveal a higher diversity of diets for copepods 
without gut dissection, which could help reconstruct copepod 
feeding conditions in natural waters. As demonstrated by previous 
studies that copepods can graze on other copepods or even their 
own nauplii to meet their metabolic costs in field, the PCR assay 
developed in this study might cause underestimation of diet 
diversity as designed to exclude copepods and thus would miss 
eggs, larvae, or adults of these copepods that might be a prey 
[2,16]. Possible retrieval of the missing information would rely on 
species-specific PCR assays requiring prior knowledge about what 
species may be a prey. For instance, Durbin et al. (2008) reported 
that adult female Centropages typicus can feed on nauplii of ^. 
tonsa by using species-specific primers for the mitochondrial 
cytochrome oxidase subunit one gene (mtCOI) [35]. 
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